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ABSTRACT 


When the GPV is under partial shading, several peaks appear in the 
characteristic P-V, namely a GMP and one or more local maximums. The 
classical algorithm ‘P&O’ MPPT cannot converge on the GMP for low 
irradiation values and is trapped by tracking down a LMP so making the 
algorithm ineffective making the algorithm ineffective, in this case under 200 
W/m?. An alternative objective function is developed to optimize the 
performance of the FLC by selecting the appropriate gains using PSO. In this 
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1. INTRODUCTION 

The use of renawble energy has become a major necessity for the future of our planet. The MPPT 
becomes complex when the photovoltaic generator is in the presence partial shading. The research was 
directed towards optimizing the performance of a solar photovoltaic system by extracting its maximum 
power for all atmospheric conditions. The GPV is formed by a series/parallel combination of photovoltaic 
modules to converge on the desired GMP. Several MPPT algorithms have been proposed by researchers [1]- 
[9], the most used are: Perturb and Observe, Incremental Conductance, Constant Voltage. These methods are 
based on finding the maximum of the carateristic P-V of the PV solar system and are effective only when the 
curve has only one maximum. When passing cloudy or obstacle on a GPV, multiple points appear on the 
power-tension characteristic with an overall GMP maximum and LMP points. Conventional methods do not 
converge on the desired GMP [10]-[13]. 


2. SIMULATION OF A GPV UNDER PARTIAL SHADING 

Partial shading has a major malfunction of a GPV under better conditions. Recent studies were 
interested in the effect of partial shading on the energy efficiency of a GPV and its impact on the functioning 
of the photovoltaic system to be studied [17], [18]. This phenomenon has pushed researchers [19]-[24] to use 
methods and algorithms to maximize power delivered by the GPV. In our previous work [25] using the 
classical 'P&O' alogrithm for rechrche of the MPPT of a photovoltaic generator under partial shading for 
different irradiation values and a constant temperature for a configuration of one hundred series/parallel 
photovoltaic modules divided into three groups (G1 to G3) each subjected to different shade levels in 
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Figure 1. We found that this approach to the 'P&O' algorithm leads to an error in tracking the global power 
point when it changes position from Pm3 to Pm2 for values under 200 W/m? in figure 2. The special 


categorization and terminology given by [26] is used for this purpose. Figure 1 shows a GPV composed of 
100 modules in a mixed group (10x10). 
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Figure 1. GPV (10x10) under parial shading divided 


Figure 2. P-V curve of GPV under partial shading of 
into three groups 








Figure 1 
The I-V characteristic of the PV module is then given by [18]. 
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3. RESULTS AND DISCUSSION 


This simulation was made for the configuration of Figure 1. The non-shaded modules receive 
irradiation level of G=lkW/m?. The shaded modules are under G=0.1kW/m?. The characteristic P-V 
illustrates the presence of three peaks with a single GMP and two local maximums in Figure 2. Figure 3 
brings together all the P-V curves for different irradiations changes in steps of 0.1kW/m? and a constant 
temperature T=25°C of the Figure 1. From the simulation of Figure 3, the power increases as the irradiation 
increases. We also note the change in the position of the GMP from Pm3 to Pm2 for low irradiation values. 


Table 1 shows the effect of partial shading on the ‘P&O’ MPPT algorithm for a grid of values of G. It gives 
the ‘P&O’ tracked power and the maximum available power. 
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Figure 3. GPV under partial shading given by Figure 1 for different shading irradiation levels 
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Figure 3 shows that for a sudden change in irradiation for shaded modules of 1000 to 50 and 100 
W/m2?, the classical "P&O" algorithm does not have the intelligence to converge the GMP. In this case the 
power losses are 08,24% and 08.99%. Figure 4 shows that for a sudden change in the shading irradiation 
level (from 1000 to 100 W/m72), the MPP tracked value using the P&O algorithm converges to 2886 W and 
not towards the global MPP who is 3208 W. The idea is to optimize the maximum power at the exit of the 
GPV, in this case GMP by inserting the FLC control whose gains are regulated by PSO in order to optimize 
the maximum output power of the GPV. 








Table 1. The ‘P&O’ algorithm Power tracking loss 
for different shading irradiation levels with the S °° “Oy 
shading configuration of Figure 1 S Tracked MPP w ith PO > 
G Maximum Power Tracked Power Loss G ais ra | 
(W/m?) (W) Power (W) (%) Sada ee tay ae |Pm3=2886 W 
50 3099 7587 08,24 E Al shaded i EN | F 
100 3208 2886 08,99 = \ Tt \ / 
200 3488 3488 0 S eae, INN 
300 4093 4093 0 5 ab Miustan = % 
400 4698 4698 0 z F a a ea Sa Y 
500 5300 5300 0 | ee a 
600 5896 5896 0 S 40 30 120 160 200 l 
700 6480 6480 0 Voltage across PV array [V] 
800 7045 7045 0 
900 7570 7570 0 . . . 
1000 7997 7997 0 Figure 4. Classical ‘P&O’ tracking process from a 


uniform standard irradiation to the partial shading 
configuration of figure 1 with a shading levelof 
0.1kW/m2 


4. SEARCH FOR THE GMP OF A GPV UNDER PARTIAL SHADING BY THE FLC 

Explaining in this section we are going to propose a FLC command to give it the ability to track the 
global MPP of a partially shaded PV array. In the presence of uniform irradiation, a single GMP appears on 
the P-V chacarteristic and it is easy to track it down by the classical 'P-O'. In the event of partial shading, the 
power output from GPV will decrease significantly and several multiple points may appear in the P-V 
characteristic. In this case, the P-O algorithm has not the potential to find the GMP unless the operating point 
does not change position. Our model is set with a limited number of rules [27]-[29]. The FLC command used 
in the search for the GMP performs better, namely its robustness and simplicity. The main parts of FLC, 
fuzzification, rule-base, inference and defuzzification, are show in Figure 5 [30]. The strategy to simulate our 
configuration of Figure 1 is to insert between the GPV and the load a DC-DC converter (boost) controlled by 
the FLC algorithm. 
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Figure 5. The FLC structure 


The principle of a FLC command is based on two input variables which are the error E and the error 
change dE and an output variable dD (variation of the duty cycle). The value of the output variable, which 
drives the static converter to find the GMP, is determined using a truth table and changing input parameters. 
error E is defined as the difference between dP/dV(k) and the desired value dP/dV (k)=0. This last value 
corresponds to the unique extremum of the curve P=f(V). This extremum is a maximum. The more positive E 
is, the higher the value of P. Conversely, the more negative E is, the lower the value of P. Finally, when E 
tends to 0, the value of P tends to its maximum, the MPP. It can be likened to the slope of P=f (V). 


E(k) ==(k) -0=— (3) 
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The variation in the deviation dE indicates in which direction and in what proportion the deviation 
changes as the algorithm proceeds. So, when dE tends to O, the system stabilizes. The first, second input 
variable and the only one output for the FLC are divided into five fuzzy sets: NB (negative big), NM 
(negative medium), Z (zero), PM (positive medium) and PB (positive big). The rules are formed as shown in 
Table 2 and Figure 7 [31]. 


dE(k) = E(k) — E(k — 1) (4) 


Table 2. Fuzzy logic rules 

dE/E NG, NM, Z, PM, PBB 
NG NG, NG, NG, NM, EZ 
NM NG, NG, NM, EZ, PM 


EZ NG, NM, EZ, PM, PG 
PM NM, EZ, PM, PG, PG 
PG EZ, PP, PG, PG, PG 














Figure 7. Input (E), (dE) and output (dD) membership Figure 8. Characteristic surface of the FLC 
functions for proposed FLC command 
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Figure 9. A fuzzy MPPT optimization approach 
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The FLC command depends heavily and necessarily on the correct selection of certain design 
parameters namely input gains and output gain to improve the performances of the structure of the FLC 
method. To guide us in searching for the appropriate gains for the FLC controller, we used the particle swarm 
optimization method to optimize the parameters of the belonging functions associated with the input and exit 
variables. The idea is to optimize our maximum power delivered by the GPV in this case GMP inserting the 
command FLC by using a PSO program for the research of the appropriate gains of the FLC applied to the 
pursuit of the MPP. PSO as one of meta-heuristic optimization techniques depends on searching the optimal 
solution within the search area based on the exchange of experiences among particles in the population. The 
particles in the swarm modify their positions in the next iteration based on individual local best position and 
global best position of swarm. Each particle represents a solution for the control variables of the entire 
optimization problem [32]. PSO is used to optimize the gains of FLC algorithm which optimize the power 
output from PV array (MPPT) [33]-[35]. The proposed system simulated in MATLAB environment. To this 
end a ‘fuzzy _optmize_PSOMATLAB' program has been designed for the search for the gains of the FLC. 
The appropriate gains for this optimization approach by running this PSO program as presented in Table 3. 


function Fuzzy Optimize PSOMATLAB 

warning off 

tic 

rng default 

options = 

optimoptions ('particleswarm', 'PlotFcns',@pswplotbestf, 'MaxIter',20,'SwarmSize',20); 
lb = 0.00005*ones (1,2); 

ub = 0.009*ones (1,2); 


[x, fval,exitflag,output] = particleswarm (@func,2,1b,ub,options); 
toc 
end 
function fitness=func (x) 
assignin('base', 'Ta', [x(1) x(2)]) 


sim('Panneau Shading MPPT',0.1); 
fitness=max (Sys Out1.Data) 
end 


Table 3. The performance of FLC 


Variable Values 


gE 0.001 
gdE 0.0001 
odV 3900 


5. COMPARATIVE SIMULATION RESULTS AND ANALYSIS 

In this FLC method the presence of partial shading is achieved by a power variation of +100W. The 
proposed solution for GMP research is to scan the speech universe [0 1] of the duty cycle by making an 
extension of the values of the FLC parameters while saving the maximum value of the power obtained in the 
process. This will allow us to detect the real GMP. A MATLAB program has been developed to simulate the 
implementation of the method proposed to the MPPT controlled GPV system in Figure 10 as part of partial 
shading. To test the FLC algorithm, the search for GMP was simulated using a GPV (10x10). The GPV is 
uniformly insolated at first, before the shading configuration of Figure 1 occurs at t=1s with a shading level 
of G=100W/m2. 

Figure 11 (a) and Figure 11 (b) shows the simulation results. When shading occurs, it appears that 

the power generated by the GPV evolves along the stabilizing PV characteristic at the GMP namely 3208 W, 
that the classical algorithm ‘P&O’ does not arrive to follow him in Figure 11 (a). The Figure 11 (b) shows us 
the ability of the FLC to track down the global power point that is 3208W that the classical 'P&O' algorithm 
can't track it. It should be noted, however, that ripples appear during the simulation by the FLC command but 
less important compared to those produced by the 'P&O' algorithm. Any attempt to delineate these ripples by 
changing the earnings values of our FLC controller leads us to move away from the global power point and 
thus to cause losses. The adjustment of the parameters (gains) of the FLC controller by the PSO algorithm 
allowed us to avoid a work of searching for gains by tweaking that can take a lot of time. 
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Figure 10. Bloc simulink of the GPV under partial shading with MPPT methods 
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Figure 11. Curve of P(t) with variation of G (1000W/m? for100W/m72) of Figure 1 with, (a) ‘P&O’ and 
(b) FLC 
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Figure 12. FLC tracking process from a uniform standard irradiation to the partial 
shading configuration of Figure 1 with a shading level of 0.1kW/m? 
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6. CONCLUSION 

In this article, we show the inability of the classical 'P&O' algorithm to track the overall power point 
for <200W/m2 values causing losses for a GVP in the presence of partial shading. We then proposed an 
optimization approach to a fuzzy logic MPPT command that moves the operating point imposed by the 
partial shading configuration to a point near the global. The results showed a perfect ability of the proposed 
algorithm to track the overall MPP for the shading configuration of Figurel for a low level of irradiation. 
Finally, these results lead us to say that the proposed FLC control method can solve the overall problem of 
tracking MPs with the classical 'P&O' algorithm. 


APPENDIX 
Table 4. PV module manufacturer’s data 
Parameters vakue 
Open circuit voltage V,, 21.9V 
Short-circuit current Isc 4.95 A 
Optimal voltage Vn 17.5 V 
Optimal current Im 4.57 A 
Maximal power Pn 80 W 
Series resistance R at T,,, 0.0102 Q 
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